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Abstract. We introduce a notion of Aubry set for weakly coupled systems of 
Hamilton-Jacobi equations on the torus and characterize it as the region where the 
obstruction to the existence of globally strict critical subsolutions concentrates. 
As in the case of a single equation, we prove the existence of critical subsolutions 
which are strict and smooth outside the Aubry set. This allows us to derive in a 
neat way a comparison result among critical sub and supersolutions with respect 
to their boundary data on the Aubry set, showing in particular that the latter is a 
uniqueness set for the critical system. Furthermore, we show that the trace of any 
critical subsolution on this set can be extended to the whole torus in such a way 
that the output is a critical solution. We also highlight some rigidity phenomena 
taking place on the Aubry set: first, the values taken by the differences of the 
components of a critical subsolution, on this set, are independent of the specific 
subsolution chosen; second, for each point y in the Aubry set, there exists a vector 
which is a reachable gradient at y of any critical subsolution. 



In the study of the Hamilton- J acobi equation for Tonelli Hamiltonians, weak 
KAM theory |14[ [15j [TBI [18] is a bridge between the PDE per se, and tools from 
the theory of dynamical systems. The use of the Lax-Oleinik formula, which is a 
variational formula to represent the viscosity solutions of the evolutionary equation, 
allows one to make rigorous the intuition that a solution of the Hamilton-Jacobi 
equation 



is deeply linked to the image of the graph of D x u(0, •) through the Hamiltonian 
flow of H. It is therefore rather natural that some invariant and minimizing sets, 
the Aubry and Mather sets, capture the long time behavior of the evolutionary 
equation. These sets are included in what can be seen as generalized invariant 
Lagrangian manifolds. They appear as graphs of differentials of weak KAM (or 
critical) solutions, that is viscosity solutions of the equation 



where c is a real number also known as critical value. 

From the PDE viewpoint the Aubry set turns out to be the region where the 
obstruction to the existence of globally strict critical subsolutions concentrates. In 
particular, the existence of a critical subsolution, smooth and strict outside the 
Aubry set O [19j [20], implies that the latter is a uniqueness set for the critical 
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equation: two critical solutions agree on the whole torus if they agree on the Aubry 
set. 

The purpose of the present paper is to generalize some PDE aspects of weak 
KAM theory to the case of weakly coupled systems of first order Hamilton-Jacobi 
equations. More precisely, we consider a family of systems of the form 

771 

Hi(x, Dm) + bij(x)uj(x) = a in T N for every i G {1, . . . , m}, (3) 

3=1 

where a is a real constant, Hi, . . . ,H m are continuous Hamiltonians defined on 
the cotangent bundle of T N , convex and coercive in the momentum variable, and 
B(x) := (bij(x)) is the coupling matrix, i.e. an m x m matrix with continuous 
coefficients satisfying 

m 

b%j{x) ^ for j 7^ i, bij(x) ^ for every x £ T N and i £ {1, . . . , m}. 

3=1 

The coupling matrix is additionally assumed irreducible, meaning, roughly speaking, 
that the coupling is non-trivial and the system cannot be split into independent 
subsystems, see Section [L2l for the precise definition; and degenerate, i.e. 

m 

bij(x) = for every x € T N and i G {1, . . . , m}. 

3=1 

Under these assumptions, there exists a unique constant a G M. for which the 
system ([3]) admits viscosity solutions. We characterize such a quantity, hereafter 
denoted by c and termed critical value, as the minimal a£l for which the corre- 
sponding weakly coupled system admits viscosity subsolutions. 

We then study the corresponding critical weakly coupled system and show that the 
obstruction to the existence of a globally strict subsolution is not spread indistinctly 
on the torus, but concentrates on a closed set A, that we call Aubry set in analogy 
to the case of a single critical equation. In particular, we show the existence of 
a critical subsolution (i.e. a viscosity subsolution of the critical weakly coupled 
system), which is strict outside the Aubry set. Furthermore, any such subsolution 
can be taken of class C°° in the complementary of A. This is achieved by exploiting 
the regularization procedure presented in [191 [20] i n the case of a single equation. As 
a byproduct, we obtain that all such subsolutions form a dense subset of the family 
of critical subsolutions. 

The analysis outlined above allows us to derive in a neat way a comparison result 
among critical sub and supersolutions satisfying suitable "boundary" conditions on 
A, see Theorem 15.31 This generalizes to our setting Theorem 3.3 in [5], therein 
established for Hamiltonians of Eikonal type. In particular, we derive that the Aubry 
set is a uniqueness set for the critical weakly coupled system. We furthermore show 
that the trace of any critical subsolution on A can be extended on the whole torus 
in such a way that the output is a critical solution, see Theorem 15.51 

Our study highlights a rigidity phenomenon taking place on the Aubry set, see 
Theorem 15. II any two critical subsolutions differ, at each point y of A, by a vector of 
the form k (1,1, ... ,1), where k is a real number depending on y and on the critical 
subsolutions itself. This means that the differences of the components of a critical 
subsolution on A are actually independent of the specific subsolution chosen. The 
proof relies on Proposition 12. 3( which is a key remark for this and other results in 

2 



the paper and which depends in a crucial way on the irreducible hypothesis assumed 
on the matrix B(x). In the particular case when there exists a critical subsolution of 
the kind (y(x),v(x), . . . , v(x)) , we infer that any other critical subsolution is of this 
form on A. This accounts for the kind of symmetry already observed, in a weaker 
form, in [5] for the particular class of Hamiltonians therein considered, see Section 
16.11 for more details. 

A second rigidity phenomenon that we point out is when the Hamiltonians are 
additionally assumed strictly convex in the momentum: in this case we prove that, 
at any point of the Aubry set, the intersection of the reachable gradients of all the 
critical subsolutions is always nonempty, see Proposition 14.41 This can be regarded 
as a weak version, for weakly coupled systems, of a result holding in the case of a 
single equation. Under suitable regularity assumptions on the Hamiltonian, it is in 
fact well known that the subsolutions of the critical equation ([2]) are all differentiate 
on the Aubry set and have the same gradient, see [181 fT9| [20] . 

We end our study by presenting a couple of situations where more explicit results 
can be obtained for the critical value and for the Aubry set, see Section [6l In the 
first example, we focus on the setting studied in [5] and we show that our notion of 
Aubry set is consistent with the one therein given. The second example contains, as 
a particular instance, the case when the Hamiltonians are all equal, say to H. In this 
specific situation, we show that the critical value and the Aubry set of the weakly 
coupled system agree with the critical value and the Aubry set of H. We furthermore 
show that the solutions of the critical system are all of the form u(x)(l, . . . , 1), where 
u is a critical solution for H, thus showing that, as far as critical solutions of the 
system are concerned, the coupling is not playing any effective role. 

Weakly coupled systems are a particular instance of monotone systems, see [U 
[TBI |2"5], 124] . but our interest for the issues herein considered stems from a series 
of recent papers [251 EH EH] addressed to understand the long-time behavior of 
the solutions (ui(t,x), . . . ,u m (t,x)) to the evolutionary counterpart of ([3]), i.e. the 
system 

du- m 

-^-+H i {x,D x u i )+^b ij {x)u j {t,x) = Q in (0, +oo) x T N , Vi G {1, . . . , m}. (4) 
3=1 

The issue is to find general conditions on the Hamiltonians and on the coupling 
matrix sufficient to guarantee that 

Ui(t,x)+ct =4 V{ in for every i 6 {1, m}, 

t— 5-+0O 

where c is the critical value and (v\(x), . . . , v m {xy\ is a critical solution. The problem 
is non-trivial due to the lack of a strong comparison principle for the critical system. 

The case of a single equation, i.e. equation i^j, is by now well understood. The 
first convergence results were obtained by Namah and Roquejoffre [27] for convex 
Hamiltonians of Eikonal type, but the breakthrough is due to Fathi |17j . maybe 
historically one of the first PDE achievements of weak KAM theory. By making use 
of the dynamical insight, Fathi established the long-time convergence for Hamilto- 
nians of Tonelli type, i.e. smooth, strictly convex and superlinear in the momentum 
variable. Related results can be also found in [29]. The dynamical approach was 
subsequently relaxed to much less regular Hamiltonians (continuous, coercive and 
strictly convex in the momentum variable) by Davini and Siconolfi [10J, then simpli- 
fied and extended by Ishii [22] to the non-compact setting. A completely different 
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approach, based on PDE methods, was instead proposed by Barles and Sougani- 
dis [2]: they establish the long-time convergence for continuous Hamiltonians that 
are coercive and satisfy a convex-type inequality with respect to the momentum 
variable. This condition includes the case of Tonelli Hamiltonians and also some 
nonconvex functions. 

Thanks to the qualitative analysis available for the critical equation, the asymp- 
totic problem can be reduced to studying the convergence on the Aubry set. The 
issue is however subtle: the afore mentioned references reveal that, in order to have 
a general convergence result, some kind of strict convexity in the momentum has to 
be required on the Hamiltonian. 

The question whether a similar statement holds for weakly coupled systems was 
first raised in [5] and answered positively by considering a setting analogous to [27] . 
More precisely, the Hamiltonians are assumed of the form 

Hi(x,p) := F{(x,p) — Vi{x) for every (x,p) € x and i € {1,... ,m}, 

where each Vi is a continuous, non-negative function on and each F{ is a con- 
tinuous function on x M^, convex and coercive in p, and satisfying 

Fi(x,p) ^ Fi(x,0) = for every (x,p) eT N x R N . 

The asymptotic convergence result is proved by assuming that the set 

m 

A := p| V; -1 ({0}) <~){xeT N : B(x) is degenerate } 

1=1 

is non-empty (note that here the coupling matrix is not assumed degenerate on the 
whole T ). In particular, this implies that the critical value c is equal to 0. The 
authors first prove that Ao is a uniqueness set for the critical weakly coupled system 
with respect to suitable "boundary" conditions on *4o- Then they exploit the special 
structure of the Hamiltonians to show that an appropriate linear combination of the 
components of the time-dependent solution asymptotically converges on Ao, and 
this is actually enough to conclude. Similar results have been also obtained in [26j . 

A significant step forward has been recently taken by Mitake and Tran [25J and 
by Nguyen [28], independently. In both papers, the coupling matrix is assumed 
independent of x and degenerate. 

In [25], the convergence is established in the case m = 2 for continuous Hamilto- 
nians that are strictly convex and coercive in the momentum variable. The analysis 
is based on an interesting dynamic programming formula for systems whose states 
are governed by random changes and on an adaptation of the techniques exploited 
in |ini[22|. 

The same result has been proved in |28] for systems with an arbitrary number 
of equations through a completely different approach based on PDE techniques. 
Furthermore, the author is able to extend the convergence result under a different 
set of assumptions in the same spirit of [2], thus including some examples of non- 
convex Hamiltonians. 

The case of a coupling matrix depending on x is widely open. The analysis carried 
out herein sheds some light on the kind of phenomena occurring for the stationary 
system at the critical level and, we hope, might be useful for further investigations 
on the subject. 

This paper is organized as follows. Section [l.ll contains the notation. In Section 
11.21 we present some basic results of linear algebra on coupling matrices. In Section 
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11.31 we precise the setting and the assumptions adopted throughout the paper, and 
we present a brief overview of existing results on weakly coupled systems. In Section 
[2] we give the definition and establish some relevant properties of the critical value. 
In particular, we prove the existence of viscosity solutions to the critical weakly 
coupled system, see Theorem 12.121 Some proofs are postponed to the Appendix [Al 
In Section [3] we give the definition of Aubry set and explore its properties. This 
is done by introducing an analogue of the Mane potential, the Mane matrix. The 
first part of Section U] is devoted to the regularization of subsolutions outside of the 
Aubry set. In the second part, we prove the rigidity phenomenon enjoyed by the 
reachable gradients of the critical subsolutions previously described. The fact that 
the values assumed on the Aubry set by a critical subsolution depend very less on the 
subsolution itself is instead proved at the beginning of Section [5j In the remainder 
of the section we state and prove the comparison principle, Theorem 15. 3\ and we 
show how the trace of a critical subsolution on the Aubry set can be extended to 
the whole torus to produce a critical solution, see Theorem 15.51 Last, in Section [6] 
we illustrate our theory on two examples. 
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Paris, that he acknowledges for the kind hospitality. The first author has been sup- 
ported by Sapienza Universita di Roma through the Research Project 2011 Analisi 
ed approssimazione di modelli differenziali nonlineari in fluidodinamica e scienza 
dei materiali. 

1. Preliminaries 

1.1. Notations. Throughout the paper, we will denote by the A^-dimensional 
flat torus, where N is an integer number. The scalar product in will be denoted 
by ( • , • ), while the symbol | • | stands for the Euclidean norm. Note that the latter 
induces a norm on T , still denoted by | • |, defined as 

\x\ := min \x + k\ for every x £ T^. 

K&Z N 

We will denote by Br(xq) and Br the closed balls in of radius R centered at xq 
and 0, respectively. 

With the symbol M+ we will refer to the set of nonnegative real numbers. We say 
that a property holds almost everywhere (a.e. for short) in a subset E of if it 
holds up to a negligible subset of E, i.e. a subset of zero iV-dimensional Lebesgue 
measure. 

By modulus we mean a nondecreasing function from M + to R+, vanishing and 
continuous at 0. A function g : M + — > R will be termed coercive if lim g(h) = +oo. 

h— >+OD 

We will say that (p n ) n is a sequence of standard mollifiers if p n (x) := n N p{nx) in 
M. N for each n € N, where p is a smooth, non-negative function on TSL , supported 
in B\ and such that its integral over M. N is equal to 1. 

Given a continuous function u on T N , we will call subtangent (respectively, su- 
pertangent) of u at xq a function <fi of class C 1 in a neighborhood U of xo such that 
4>(xo) = u{xq) and <f>{x) ^ u(x) for every x € U (resp., ^). Its gradient D<f>{xo) 
will be called a subdifferential (resp. superdifferentiat) of u at xq, respectively. The 
set of sub and superdifferentials of u at xq will be denoted D~u(xq) and D + u(xo), 
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respectively. We recall that u is differentiable at xq if and only if D + u{xq) and 
D~u{xq) are both nonempty. In this instance, D + u{xq) = D~u(xq) = {Du(xq)}, 
where Du{xq) denotes the differential of u at xq. We refer the reader to [6] for the 
proofs. 

When u is locally Lipschitz in T N , we will denote by 8*u(xq) the set of reachable 
gradients of u at xq, that is the set 

d*u(xo) = {lim Du(x n ) : u is differentiable at x n , x n — > xo }, 

n 

while the Clarke's generalized gradient d c u(xo) is the closed convex hull of 8*u(xq). 
The set d c u(xo) contains both D + u{xq) and D~u(xq), in particular Du(xq) S 
d c u(xo) at any differentiability point xq of u. We refer the reader to [8] for a detailed 
treatment of the subject. 

We will denote by ||<?||oo the usual L°°-norm of g, where the latter is a measurable 
real function defined on T^. We will write g n g in to mean that the sequence 
of functions (g n )n uniformly converges to g in T N . We will denote by (C(T Ar )) the 
Banach space of continuous functions u = (u\, . . . ,u m ) T from to M m , endowed 
with the norm 

||u-v||oo= max \\ui-Vi\loo, u, v € (C(T 7V )) m . 

We will write u n u in to mean that 1 1 u."' — u||oo — > 0. A function u € (C(T iV )) m 
will be termed Lipschitz continuous if each of its components is K-Lipschitz contin- 
uous, for some k > 0. Such a constant k will be called a Lipschitz constant for u. 
The space of all such functions will be denoted by (Lip(T Ar )) . 

We will denote by 1 = (1, • • • , 1) T the vector of M m having all components equal 
to 1, where the upper-script symbol T stands for the transpose. We consider the 
following partial relations between elements a, b G R m : a ^ b if a\ ^ b{ (resp., <) 
for every i € {1, . . . , m}. Given two functions u, v : — > R m , we will write u ^ v 
in (respectively, <) to mean that u(x) ^ v(x) (resp., u(x) < v(x)) for every 
x G T N . 

1.2. Linear algebra. Here we briefly present some elementary linear algebraic re- 
sults concerning coupling matrices. 

Definition 1.1. Let B = (bij)ij be a m x m-matrix. 

(i) We say that B is a coupling matrix if it satisfies the following conditions: 

m 

bij < for j ^i, ^2 b{j > for any i 6 {1, ... , m}. (C) 
3=1 

Lt is additionally termed degenerate if 

m 

bij = for any i = 1, . . . , m. 

3=1 

(ii) We say that B is irreducible if for every subset X C {1, . . . ,m} there exist 
i Gl and j Z such that bij ^ 0. 



When a coupling matrix is also irreducible, further information can be derived on 
its elements. We have 
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Proposition 1.2. Let B = (bij)ij be an irreducible m x m coupling matrix. Then 
bu > for every i € {1, . . . ,m}. 

Proof. Indeed, if bi i = for some io £ {1> ■ • ■ > m }> condition ((C|) would imply 
bi j = for every j £ {1, . . . , m}, in contradiction with the fact that B is irreducible. 

□ 

The following invertibility criterion holds: 

Proposition 1.3. Let B = j be an m x m irreducible coupling matrix. Then 

(i) Ker(B) C S pan{(l, . . . , 1) T } = Ml; 

(ii) Ker(i?) = span{(l, . . . , 1) T } = Ml if and only if B is degenerate. 
In particular, B is invertible if and only if 

m 

bij > for some i £ {1, . . . , m}. 

3=1 

Proof. We first remark that, by assumption (fU|) . 

bu ^ ^2 \bij | for every i G {1, . . . , m}. (5) 

Let us prove (i). Let v = (v\, . . . ,v m ) T £ Ker(B) and set 

I={iG{l,...,m} : Vi = m&x{vi, . . . ,v m } }. 

We claim that X = {1, . . . , m}. Indeed, if this were not the case, by the irreducible 
character of B there would exist i £ X and X such that ^ 0. Since £>v = 0, 
we would get in particular 

buvi = ^2vj\bij\ < «j ^ |6y| < Uj6ji. 

Then the inequalities must be equalities. We infer 

Vj\bij\ = Vi\bij\ for every j / i, 

in particular Vk = Vi = max{ui, . . . , v m }, yielding that k belongs to X, a contradic- 
tion. 

The remainder of the statement trivially follows from item (i). □ 

The following proposition gives an obstruction to being in the image of a degen- 
erate coupling matrix. 

Proposition 1.4. Let B a coupling and degenerate mxm matrix. If a = (ai, . . . , a m ) 
satisfies > for every i 6 {1, . . . , m}, then a ^ lm(B). 

Proof. Let us assume by contradiction that there exists v = (vi, . . . ,v m ) T such 
that 

Bv = a. 

Let Vk = min{t>i, . . . ,v m }. We have 

m m 

a k = ^2 bk 3 v i ^ hk i Vk = °' 

3=1 3=1 

in contradiction with the hypothesis at > 0. □ 
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1.3. Weakly coupled systems. Throughout the paper, we will call convex Hamil- 
tonian a function H satisfying the following set of assuptions: 

(HI) H : T N x R N -> R is continuous; 

(H2) p t — y H(x,p) is convex on R^ for any x G T^; 

(H3) there exist two coercive functions a, f3 : R + — > R such that 

a (\p\) < H(x,p) ^ (\p\) for all (x,p) £T N x R N . 

The Hamiltonian H will be termed strictly convex if it additionally satisfies the 
following stronger assumption: 

(H2)' p i — y H(x,p) is strictly convex on R^ for any x G T . 

Moreover, we will denote by B(x) = (bij(x)}.. an m x m-matrix with continuous 

coefficients bij(x) on T N . If not otherwise stated, the following hypotheses will be 
always assumed: 

(Bl) B{x) is an irreducible coupling matrix for every x G T^; 

(B2) B{x) is degenerate for every x G T N . 

Let Hi(x,p), . . . ,H m (x,p) be convex Hamiltonians, i.e. functions satisfying con- 
ditions (H1)-(H3). We are interested in weakly coupled systems of the form 

Hi(x, Dui) + (.B(x)u(x)) i = dj in for every i E {1, . . . ,m}, (6) 

for some constant vector a = (ai, . . . , a m ) T , where u(x) = (ux(x), . . . , n m (x)) T and 
(S(x)u(a;)) . denotes the i-th component of the vector £>(x)u(x), i.e. 

m 

(B(x)u(x)) i = ^2b ij (x)u j (x). 
i=i 

Remark 1.5. The weakly coupled system flSJ) is a particular type of monotone 
system, i.e. a system of the form 

Gi(x,ui(x), . . . ,u m (x), Dui) = in for every i G {1, m}, 

where suitable monotonicity conditions with respect to the tij-variables are assumed 
on the functions Gi, see [U HU Ell ESJ El]. In the specific case considered in this 
paper, the conditions assumed on the coupling matrix imply, in particular, that each 
function Gi is strictly increasing in Uj and non-increasing in Uj for every j ^ i. This 
kind of monotonicity will be exploited in many points of the paper. 

Let u e (C(T 7V )) m . We will say that u is a viscosity subsolution of © if the 
following inequality holds for every (x,i) eT N x {l,...,m} 

Hi(x,p) + (-B(x)u(a;)) . ^ a, for every p € D + Ui{x). 

We will say that u is a viscosity supersolution of (J6]) if the following inequality holds 
for every (x, i) G x {1, . . . , m} 

Hi(x,p) + (f?(x)u(x)) . ^ ai for every p G D~Ui(x). 

We will say that u is a viscosity solution if it is both a sub and a supersolution. In 
the sequel, solutions, subsolutions and supersolutions will be always meant in the 
viscosity sense, hence the adjective viscosity will be omitted. 
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Due to the convexity of the Hamitonian Hi, the following equivalences hold: 

Proposition 1.6. Let a G R, i G {1, . . . ,m} and u G (Lip(T Ar )) m . 27ie following 
facts are equivalent: 

(i) Hi(x,p) + [B(x)u(x)^ i ^ a for every p G D + U{(x) and x G T^; 

(mJ Hi(x,p) + (S(aj)u(aj)) . ^ a for every p G D~Ui(x) and x G T^; 

(mj Hi(x,p) + (S(aj)u(aj)) . ^ a for every p G d c Ui{x) and x G T^; 

(wj Hi(x, Duiixfj + (B(aj)u(aj)) . ^ a /or a.e. x G T^. 

Next, we state a proposition that will be needed in the sequel, see also [131 122 
[2T| [23] for similar results. 

Proposition 1.7. Lei F be a subset of {C{T N )'\ m and define the functions u, u on 
by setting: 

Ui(x) = inf Ui(x), Ui(x) = sup U{(x) for every x G and i G {1, . . . , m}. 

Assume that u and u belong to (C(T Ar )) m and Zei a G M. m . Then: 

(i) if every u £ J 7 is a subsolution of ©, then u is a subsolution of ©; 

(ii) if every u G T is a supersolution of ©, then u is a supersolution of ((6]). 

The previous proposition is analogous to a well known fact for scalar Hamilton- 
Jacobi equations, see for instance Section 2.6 in pQ. The proof can be easily recovered 
by arguing similarly and by exploiting the monotonicity of the system. 

We will be also interested in the evolutionary counterpart of ([6]), i.e. the system 

-^- + H i (x,D x u i ) + (B(x)u(t,x)) i = in (0, +oo) x V* G {1, . . . , m}, (7) 

where we have denoted by u(t, x) = (ui(t, x), . . . , u m (t, x)^ T . 

The following comparison result holds, see for instance [1] for a proof. 

Proposition 1.8. Let T > and v, u G (Lip([0, T] x T N )^ m be, respectively, a sub 
and a supersolution of (J7]). Then, for every i G {1, . . . , m}, 

Vi(t,x) — uAt,x) ^ max max (WO, •) — «i(0, •)) , (t,x) £ [0,T] x T N . 

By making use of this proposition and of Perron's method, it is then easy to prove 
the following 

Proposition 1.9. Let uo G (Lip(T Ar )) m . Then there exists a unique function 
u(t,x) in (Lip(R + xT^))™ that solves the system (|7|) subject to the initial condition 
u(0, x) = uo(x) in T N . Moreover, the Lipschitz constant ofu(t,x) in M + x only 
depends on the Hamiltonians H\, . . . ,H m and on the Lipschitz constant of uq. 

We will denote by S(t)uo(x) the solution u(i, x) of (JT|) with initial datum uo. 
This defines, for every t > 0, a map 

S(t) : (Lip(T")) m -»■ (Up(T N )) m . 

We summarize in the next proposition the properties enjoyed by such maps, which 
come as an easy application of the above results. 
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Proposition 1.10. For every t, s > and u, v € (Lip(T Ar )) m we have: 

(i) (Semigroup property) S(s)(S(t)u) =S(t + s)u inT N ; 

(ii) (Monotonicity) if v ^ u in T N , then S(t)v < S(t) u in T : 
(ni) (Non expansiveness property) ||<S(i)v — S^uHoo ^ |[v — u|joo; 
(iv) for every aeR, S(t)(u + ol) = S(t)u + at in T N . 

The fact that the coupling matrix B(x) is everywhere degenerate is crucial for 
assertion (iv). 

2. The critical value 

The purpose of this section is to define the notion of critical value for weakly 
coupled systems and to prove some relevant properties of the corresponding critical 
system. 

We start by proving some a priori estimates for the subsolutions of a weakly 
coupled system of the form ©. The following notation will be assumed throughout 
the section: 

IH = mm Hi(x,p) for each i € {1, . . . ,m}, (J. = mm //j. 

(x,p) iG{l,...,m} 

Proposition 2.1. Let a = (ax, ... , a m ) T € R m and u € (C(T Ar )) m such that 

(i3(x)u(x)). ^ dj for every x G T N and i £ {I, . . . ,m}. (8) 

Then there exists a constant M a only depending on a and B(x) such that 

(i) \\ui - UjWoa < M a for every i, j G {1, ... ,m}; 

(ii) | (-B(x)u(x)) .| ^ M a for every x £ T N and i £ {1, ... ,m}. 
Proof. It suffices to prove the assertion for a = a 1 . Let us set 

= mm mm ba(x), B* = max max 

Such quantities are finite valued. Moreover, B± is strictly positive in view of Propo- 
sition 11.21 and of the fact that B(x) is, for every x S T^, an irreducible coupling 
matrix with continuous coefficients. 

Let us now fix x S and assume, without any loss of generality, 

u\(x) < u 2 (x) < ■ ■ • ^ u m (x). (9) 

m 

First notice that, by subtracting Yl b m j(x)u m (x) = from both sides of equation 
j8|) with i = m, one gets 



-b mj (x)(u m (x) - Uj(x)) ^ a, 



yielding 



Um{x) — maxiij(x) ) > —b m j(x) ^ a. 



By exploiting Q and the degenerate character of the matrix B(x) we get 

^ u m (x) - u m -i(x) ^ - ' r < — . (10) 
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This proves assertion (i) when m = 2. To prove it in the general case, we argue by 
induction: we assume the result true for m and we prove it for m + 1. To this aim, 
we restate equation (JSj) as 

m— 1 

^ bij(x)uj(x) + (b im (x) + b im+1 (xfju m (x) + 6 im+ i(x) (u m+1 (x) - u m (x)j ^ a, 

then we exploit (TlOj) to get 

^ bij(x)uj(x) + (&i m (;c) + 6 im+ i(x)Ju m (a;) < a |l + ^-J (11) 

for every j € {1, . . . ,m + 1}. The irreducible character of B{x) applied to the set 
X = {m, m + 1} implies that 

hm(x) + b im+ i(x) > 

for either % = m or i = m + 1, let us say i = m for definitiveness. Assertion (%) 
now follows by applying the induction hypothesis to the system given by (llip with i 
varying in {1, . . . , m}, the corresponding coupling matrix being still irreducible and 
degenerate. 

To prove (ii) it suffices to note that, for every i £ {1, . . . , m}, 
- (B(x)u(x)) i = -bu{x)ui(x) + } j ( - bij(x))uj(x) 

^ -bu(x)ui(x) + ^ -6jj(x)(ui(x) + |[u; - itjUoo) 

^ (m- 1) f3*\\ui - UjWoo, 

and the assertion follows from (i) and from hypothesis ([8]). □ 

As a consequence, we derive the following result: 

Proposition 2.2. Let u = (ui, . . . ,u m ) T € (C(T Ar )) m 6e a subsolution of Q /or 
some a G M m . T/ien i/iere exists costants C a and K a , onZy depending on a, and 
i/ie coupling matrix B(x), such that 

(i) \\ui - uj ||oo ^ C a for every i, j £ {1, ... ,m}; 

(ii) u is K a -Lipschitz continuous in T . 

Proof. For each i € {l,...,m}, the following inequalities hold in the viscosity 
sense: 

(i + (B(x)u(x)) i ^ Hi{x, Dm) + (B{x)u(x)). ^ a { in 

yielding 

(B(x)u(x)) . ^ aj — /i for every x € T N . 
In view of Proposition 12.11 we get (i) and 

| (i?(x)u(x)) I ^ M a for every x € T^. 

Plugging this inequality in ([6]) we derive that Ui is a viscosity subsolution of 

Hi(x, Dui) sC 0i + M a in 

and assertion (ii) follows as well via a standard argument that exploits the coercivity 
of Hi(x,p) in p. □ 
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Next, we establish a remarkable property of weakly coupled systems. 

Proposition 2.3. Assume that v, u € (C(T Ar )) m are, respectively a sub and a 
supersolution of the weakly coupled system ([6]) for some a G M. m . Let xq G be 
such that 

v i( x o) ~ Ui(xo) = M := max max (vi — Ui) for some i G {1, . . . , m}. 
Then v(xo) = u(xo) + Ml. 

Proof. In view of Proposition 12.2^ we know that v is Lipschitz continuous. Set 

X= {ie{l,...,m} : (vi(x ) - Ui(x )) = M }. 

We want to prove that I = {1, . . . ,m}. Indeed, if this were not the case, by the 
irreducible character of the matrix B{xq) there would exist i el and k ^LT such 
that 

hk(x ) < 0. 

We now make use of the method of doubling the variables to reach a contradiction. 
For every e > 0, we set 

ip £ (x, y) = Vi(x) - Ui(y) ^ j^—, x, y G . 

Let M £ = max ip e and denote by (x £ , y £ ) a point in T N x T N where such a maximum 

X AT xT JV 

is achieved. By a standard argument in the theory of viscosity solution, see for 
instance Lemma 2.3 in pQ, the following properties hold: 

x £ ,y £ — > xq, — — — > 0, \x £ — xq\ — >0 as e — > 0. (12) 



Furthermore, 

p' £ := £ g2 6 G D~Ui(y £ ), p £ := p' £ - (x £ - xo) G D + Vi(x £ ) for every e > 0. 

By the Lipschitz character of V{ we derive that the vectors { p £ : e > } are equi- 
bounded, hence, up to subsequences and in view of the estimates (fT2|) . we infer 

p E , p' £ ^Po as e ->• 

for some vector po € K^. We now use the fact that v and u are a sub and superso- 
lution of (fT6j) . respectively, to get 

Hi{x £ ,p E ) + (B(x e )v(a: 6 )) i s$ 0, 

if,(y £ ,p;) + (i?(y e )u(y £ )). ^ 0. 

By subtracting the above inequalities and by passing to the limit for e — > we end 
up with 



<0, 



B(xq)(v(xq) - u(x ))j 
that is, since i 6l and the matrix B(xq) is degenerate, 

Mbu(xo) ^ ^2 \bij(xo)\(yj(x ) - Uj(x )) < M ^ |by(«o)| = Mbu(xo). 

Hence the above inequalities are equalities, in particular v^Xq) — Ufc(xo) = M since 
bik(xo) 7^ 0, in contrast with the fact that k g" 1. □ 
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Definition 2.4. For every a G R m , we denote by %(a) the set of subsolutions of 
the weakly coupled system Q. We will more simply write l~L(a) whenever a = al 
for some constant a € K, . 

Lemma 2.5. The sets 7i(a) are convex, closed (with respect to the topology of 
uniform convergence) and increasing with respect to the partial ordering on M m . 

Proof. Convexity and monotonicity are straightforward. The fact that the H(b) 
are closed is a direct consequence of stability of viscosity subsolutions. □ 

We now focus our attention to the case a = al. As a direct consequence of the 
definition of the semigroup S(t), we get the following assertion: 

Proposition 2.6. Let a£f and u G (LipCT^)) . Then u is a viscosity solution of 
([6]) with a = a 1 if and only if 

u = S(t)u + tat in T N for every t > 0. 
We have the following characterization: 

Proposition 2.7. Let a£l. The following facts are equivalent: 

(i) u G H(a); 

(ii) the map t h-> S(t)u + tat is non-decreasing on [0, -fee). 

In particular, the sets H(a) are stable under the action of the semigroup S(t), in the 
sense that S(t)(H{a)) C H(a). 

The proof of this proposition is rather technical and it is postponed to the Ap- 
pendix [Aj 

Definition 2.8. The critical value c of the weakly coupled system ([6]) is defined as 

c = inf{a G K : H(a) ^ }. (13) 

The following holds: 

Proposition 2.9. The critical value c is finite and H(c) ^ 0. 

Proof. By the growth assumptions on the Hamiltonians Hi it is easily seen that 
the function u = (0, . . . , 0) T is a subsolution of ([6]) for a$t with ao G R big enough. 

Let us proceed to show that c is finite valued and that Ti{c) ^ 0. Let (a n ) n be 
a decreasing sequence converging to c and let u n G H(c n ) for each n G N. Up to 
neglecting the first terms, we can assume that a n ^ ao for every n G N. Arguing as 
in the proof of Proposition ^. 2\ we obtain that the following inequalities are satisfied 
in the viscosity sense: 

Hi ^Hi{x, Duf) ^a n + M ao in 

for every i G {l,...,m} and n G N, showing that c is finite. We now exploit 
Proposition 12.21 by the monotonicity of the sets T~L(a) with respect to a, we infer 
that the functions u ra are equi-Lipschitz. Up to subtracting a vector of the form 
k n t to each u ra , we can furthermore assume that u"(0) = for every n G N, yielding 
sup n \\ui | |oo ^ L for some L G K by the equi-Lipschitz character of the sequence. 
Moreover, 

||n" — «™ ||oo ?s C ao for every j G {1, . . . , m} and n G N, 

yielding 

||w™||oo ^ C aQ + L for every j G {1, . . . , m} and n G N. 
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Up to subsequences, by the Arzela-Ascoli theorem, we infer that 

u" =J u in 

and u G H(c) by stability of the notion of viscosity subsolution. □ 

We now proceed to show that the stationary equation © possesses solutions of 
the corresponding weakly coupled system if and only if a equals the critical value c. 
We start with a preliminary result. 

Proposition 2.10. Let B(x) be a continuous irreducible coupling matrix on T N and 
let us assume that B(x) is invertible for every x G T^. Let v, u € (C(T Ar )) m be, 
respectively, a sub and a supersolution of the weakly coupled system (|6|), for some 
aeR m . Then 

v(x) ^ u(x) for every x G T . 

Proof. Arguing as in the proof of Proposition 12. 2\ we easily see that v is Lipschitz. 
Set 

M = max max (v{ — uA 

We want to prove that M ^ 0. Assume by contradiction that M > and pick a 
point xq G where such a maximum is attained. Set 

I={i6{l,...,m} : (vi(x ) - Ui(x )) = M }. 

Arguing as in the proof of Proposition 12.31 we infer that 

B(xo){v(xq) — u(a?o))^ ^ for every i G X. (14) 

If X = {l,...,m}, inequality (|14p is in contradiction with the fact that B(xq) is 
invertible, in view of Proposition 11.31 If X ^ {1, . . . , m}, we choose i £ X and k ^X 
such that bik{xo) < 0. From (fl"4"j) and the assumption that M > we infer that 

Mbii(xo) ^ ^2 \bij(xo)\(vj(x ) - Uj(x )) ^ M ^ [6^ (^n )] ^ M6„(aj ), 

which implies that Vk(xo) — Uk(xo) = M, in contrast with the fact that k £ X. □ 

Next, we prove that solutions to a weakly coupled system of the kind ([6]) with 
a = al may exist only if a equals the critical value. 

Proposition 2.11. Let a, b G M. and v, u G (C(T N )^ m such that the following 
inequalities are satisfied in the viscosity sense: 

Hi(x,Dvi) + (B(x)v(x)). < a in T N 

Hi(x,D Ui ) + (B(x)u(x)). ^ b inT N 

for every i G {1, . . . , m}. TTien b ^ a. 

Proof. Let us assume by contradiction that b > a. Up to replacing v with v + kl 
with k > big enough, we can assume 

v > u in T N . 

Let e > such that b — e > a + e. By continuity of the functions v and u, we can 
find A > such that 

||A Vi[|oo, IIAujlloo < e for every i G {1, . . . ,m}. 
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Then the following inequalities hold in the viscosity sense in T . 

H i {x,Du i ) + ((B{x) + Xl)u{x)). > b-e > a + e > H t (x, D Vi ) + ((B(x) + AI)v(x)) r 

For ever x G T , the matrix B(x) + A I is irreducible, satisfies (C) and the sum of 
the elements of each of its rows is strictly positive, hence it is invertible in view of 
Proposition 11.31 By Proposition 12.101 we conclude that 

v^u in T , 

achieving a contradiction. □ 

We are now able to prove a weak KAM theorem, following the lines of Fathi [18J. 
This result has been already obtained in literature in similar settings by making use 
of the so called ergodic approximation, see [26], [5] . 

Theorem 2.12. There exists a function u G H(c) that solves the weakly coupled 
system 

Hi(x, Dui) + [B(x)u(x)) i = c inT N for every i G {1, ... ,m} (15) 
in the viscosity sense. 

Proof. We have already proved in Proposition 12.91 that %{c) ^ 0. Let us introduce 
the quotient space H = 7i(c)\Ml. Arguing as in the proof of Proposition 12.9} it is 
easily seen that Tl is compact for the topology of uniform convergence. Indeed, it 
is isomorphic to the subset of T-L(c) of subsolutions whose first component vanishes 
at the point x = 0. Moreover, since the viscosity semigroup commutes with the 
addition of vectors of the form Al and leaves 7i(c) stable, it induces a continuous 
semigroup, denoted S, on %. 

By the Schauder-Tychonoff fixed point theorem (see [32]), S possesses a fixed 
point, that is, there exists an element u G H such that 

Vt ^ 0, S(t)u = u. 

Lifting these relations to H(c), we get 

Vt ^ there exists c% G M. such that S(t)u = u + q1, 

where u is any element in the equivalence class of u. Since S is a semigroup, one 
readily realizes that the following relations are verified: 

q +s = a + c s for every t, s > 0. 

Since S is continuous, we necessarily deduce that q = — tc for all t > for some 
constant c G R. 

The identity S(t)u = u — tcl, for all t ^ 0, implies that u is a viscosity solution of 
(I15D with c in place of c, see Proposition 12.61 But then c = c in view of Proposition 
I2.11l and the statement is proved. □ 

3. The Aubry set 

In this section we start our qualitative analysis on the critical weakly coupled 
system 

Hi(x, Dui) + (B(x)u(x)) i = in for every i G {1, . . . , m}. (16) 

Here and in the remainder of the paper we assume that the critical value c defined 
via (]13p is equal to 0. This renormalization is always possible by replacing each 
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Hi with Hi — c. Solutions, subsolutions and supersolutions of the weakly coupled 
system (|16p will be termed critical in the sequel. The family of critical subsolutions, 
we recall, is denoted by H(0). 

Our qualitative analysis on the critical weakly coupled system is based on the 
notion of Mane matrix, defined in analogy with that of the Mane potential. 

Definition 3.1. For all (x, y, i,j) € T N x T N x {1, . . . , to} x {1, . . . , m}, we define 

®i,j(x,y) = sup Vi(y)-Vj(x). 
veH(o) 

The following properties hold: 

Proposition 3.2. The Mane matrix verifies the following properties: 

(i) it is everywhere finite and Lipschtiz continuous; 

(ii) $. t j(y,-) € TL(0) for every (y,j)eT N x {l,...,m}; 
(Hi) for every (y,j) € T N x {1, . . . , m} and v € rl(0), 

v-Vj(y)l < $.j(y,0 inT N , 

namely $.j(y, •) is the maximal critical subsolution whose j-th component 
vanishes at y; 

(iv) the entries of the Mane matrix are linked by the following triangular inequal- 
ity: 

<&i, k {x, z) < *j,fc(sc, y) + ®i,j(y, z) 
for every i,j, k € {1, . . . , to} and x,y,z E T . 

Proof. The fact that the Mane matrix is well defined directly follows from Propo- 
sition 12.21 Lipschitz continuity comes from the equi-Lipschitz character of critical 
subsolutions. 

The second assertion comes from the fact that •) is, for every fixed (J, y), a 

supremum of critical subsolutions, hence itself a critical subsolution by Proposition 

o 

The third point is a direct consequence of the definition. 

The last point comes from the fact that <&.j(y, •) is the greatest subsolution whose 
j-th component vanishes at y. Since $.,fe(x, •) — $>j t k(x,y)l is a subsolution whose 
j-th component vanishes at y we obtain that 

$., fc (av) - $ jtk (x,y)l < *.j(j/,-) 5 
which is the triangular inequality to be proved. □ 

As in the case of a single critical equation, the Mane vectors are "almost" critical 
solutions, in the sense precised below: 

Proposition 3.3. Let y$ £ and iq G {1,...,to}. Then the function u = 
$.,i (yO) •) satisfies 

H io (x,Du io ) + (B(x)u(x)) iQ =0 to T N \ {y }. 

and 

Hi(x,Dui) + (B(x)u(x)) . = in T N for every i ^ io 
in the viscosity sense. 
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Proof. We argue by contradiction, following the classical argument of [18] for the 
classical Mane potential. 

Let (i,y) be such that either i ^ iq or y ^ y$. Let us assume that the viscosity 
supersolution condition is violated at (i, y). This means that there exists a C 1 
function ip such that ip{x) ^ <&ii (yo,x) for all x, with equality if and only if x = y, 
and 

Hi(x,Dil>(y)) + (B(y)$. Ao (y ,y)). < 0. 

Since ip is C , and B(-) and ^.^(xq,-) are continuous, it is clear that this strict 
inequality continues to hold in a neighborhood of y. We infer that it is possible to 
find e > small enough such that the function W{ := max{$j j j (yo ) O>'0 + £ } verifies 

Hi(x,D Wi (x)) + (B(x)w(x)). < for a.e. x € 

where w is the vector whose i-th coordinate is W{ and whose other coordinates are 
those of $. t i (yo, ■). In the case when i = iq and y ^ yo, we choose e > small 
enough in such a way that Wi(yo) = &i,i Q (yo,yo) = 0. Moreover, for every j ^ i, 

Hj(x,D Wj (x)) + (B{x)w(x)). ^ for a.e. x G T N , 

as it is easily seen from the fact that bji(-) ^ in and W{ ^ ®i,i {yoi ■)• 

We have thus shown that w is a critical subsolution with Wi Q (yo) = 0, w ^ 
$i,io(2/o,0 and w # *i,io(2/o,0) thus contradicting the maximality of $i,i Q (yo,-) 
amongst subsolutions whose io-th coordinate vanishes at yo- D 

Next, we show a strong invariance property enjoyed by the rows of the Mane 
matrix. 

Proposition 3.4. Let i,j G {1, . . . , m} and y G T . //"$. •) is a critical solution 
on T N , then $>.j(y, •) is too. 

Proof. Let us set v := &.,j(y, •) and u := <&.,i(y, •) + &ij(y, y) 1. In view of 
Proposition 13.31 we only need to show that 

H j(y,p) + { B {y)^{y))j > for every p G D~Vj(y). 

According to Proposition [3721 v ^ u in and = U{(y). The functions v and u 
being respectively a critical subsolution and a solution, we can apply Proposition ^. 31 
to infer that v(y) = u(y). This also implies that D~Vj(y) C D~Uj{y). Exploiting 
again the fact that u is a critical solution we finally get 

< Hj(y,p) + (B(y)u(y)). = Hj(y,p) + (B(y)v(y)). for every p G D~Vj(y). 

□ 

In view of the previous proposition, the following definition is well posed: 

Definition 3.5. T7ie Aubry set A for the weakly coupled system (|16p is i/ie se£ 
defined as 

A = {y G T w : ^.^(y, •) is a critical solution} , 
where i is any fixed index in {1, ... , m}. 

By the continuity of the Mane matrix and the stability of the notion of viscosity 
solution, it is easily seen that A is closed. The analysis we are about to present will 
show that the Aubry set is nonempty: as in the corresponding critical scalar case, 
we will see that A is the set where the obstruction to the existence of globally strict 
critical subsolutions concentrates. 
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Definition 3.6. Let v G H(0). We will say that is strict at y £ i/ t/iere exist 
an open neighborhood V of y and 5 > sitc/i t/iat 

Hi(x, Dvi(x)) + (i?(x)v(x)) i < —5 /or a.e. x € V. 

VFe wiZZ say fZiai is strict in an open subset U of T N if it is strict at y for every 
yeU. 

We start by establishing an auxiliary result that will be needed in the sequel. 

Lemma 3.7. Let w G "H(O) such that Wi is strict at y G T N . Then there exists 
w G %(0) such that Wi is C°° and strict in a neighborhood of y. 

Proof. By hypothesis, there exist r > and 5 > such that 

Hi(x, Dwi(x)) + (i?(x)w(x)) i < —5 for a.e. x £ £?2r(y)- 

Let </> : — >• [0, 1] be a C°°-function, compactly supported in B r {y) and such that 
</> = 1 in B r / 2 (y)- Let us denote by k a Lispchitz constant for the critical subsolutions 
and by ui a continuity modulus of i^j in x Br for some fixed R > k + || -D^Hoo- 
Let (/9 n ) n be a sequence of standard mollifiers on 1^ and define 

1p n (x) = (pn * Wi)(x) + \\p n * Wi - Wt||oo, Z G T^. 

Note that ^ Wi in for every n £ N and 

:= HV'n — ^illoo — >• as n — > +oo. 

Up to neglecting the first terms, we furthermore assume that all the d n are less than 
1. For every n G N, we define a function w n G (Lip(T Ar )) by setting 

Wj(x)=Wj(x) ifj/i, (a;) = (f>(x)ip n (x) + (l - 0(x))it;i(x) 

for every x G . It is apparent by the definition that wf is of class C°° in B r j2 (y) ■ 
Moreover the functions (wf ) n , and hence the (w n ) n , are equi-Lipschitz. Indeed, for 
almost every x G T , 

Dwf(x) = <f>{x)D^ n {x) + (1 - <j>{x))Dwi{x) + ty n (x) - Wi (x)) D</>(x) (17) 

that is 1 1 | |oo ^ k + H-D^Hoo. We want to show that n can be chosen sufficiently 
large in such a way that w n G "H(0) and 

Hi(x,Dwf(x)) + (£(x)w n (x)). < -~6" for a.e. x G S r (y). (18) 

We first note that, since wf ^ Wi and 6 3 j ^ in for every j ^ i, we have 

Hj(x, Dw](x)) + (B(x)w n (x)) j ^ in for every j ^ i. (19) 

Moreover, since w n agrees with w outside B r (y), in order to show that w n satisfies 
(fT9|) also for j = i, it will be enough, by the convexity of Hi, to prove (fT8j) . 
To this aim, we start by noticing that 

Hi(x,Dwf(x)) ^<f>(x)Hi(x,Lhf> n (x)) + (l - <f>(x))Hi(x,D Wi (x)) +u(d n \\D<l>\\ 0O ) 

(20) 
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for almost every x G T N , in view of (|17p and of the convexity of Hi. By Jensen's 
inequality, for every n > 1/r and every x G B r we have 



Hi(x,Dij) n (x)) = Hi\x, Dwi(x - y)p n (y)dy 



Hi(x,Dwi(x - y))p n {y) dy 

Bl/n 



< w(l/n)+/ Hi(x -y,Dwi(x -y))p n (y)dy 

JBy n 

< -/ - 2/)w(z - y))iPn{y)dy - 5 + w(l/n) 

JBi/„ 

< -(B(x)w n (x)) i -5 + o;(l/n)+£ n , (21) 



where 



e n := sup ||(B(. + *) w (- + z) L- 

|z|^l/n 

Since w n z| win and all these functions are equi-Lipschitz, it is easily seen that 
lim n e n = 0. Furthermore 

Hi(x,Dwi(x)) < -(£(a;)w n (x)). -<5 + e n for a.e. s G B r (y). (22) 

We now choose n > 1/r sufficiently large such that 

w (dn H-D0IU) + w(l/n) + e„ < - 

and plug (j2Tj) and ([22]) into (f20l) to finally get (fT8l) . The assertion follows by setting 
w := w n for such an index n. □ 

The next proposition shows that the i— th component of any critical subsolution 
fulfills the supersolution test on A. 

Proposition 3.8. Let y G A. Then, for every i G {1, . . . , m} and w G TL(0), 

Hi(y,p) + (-B(y)w(y)) . = /or every p G D~Wi{y). (23) 

Proof. Pick w G %(0) and set u = $>. t i(y, •) + tUj(y)l. According to Proposition 
I3.2} w ^ u and, by definition of u, Wi(y) = Ui(y), in particular D~Wi(y) C D~m(y). 
Now we exploit the fact that u and w are a critical solution and subsolution, re- 
spectively: from Proposition 12.31 we infer that w(y) = u(y), while Proposition 11.61 
implies 

> Hi(y,p) + (B(y)w(y)) . = fl^p) + (S(y)u(y)) . ^ Vp G 
Hence all the inequalities must be equalities and the statement follows. □ 
A converse of this result is given by the following 

Proposition 3.9. Let i G {1, . . . ,m}. The following facts are equivalent: 

(i) y£A; 

(ii) there exists w G Ti(0) such that is strict at y. 
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Proof. Let us assume (i). Since y G" A, the supersolution test for $. t i(y, •) is violated 
at (i,y). This means that there exists a C 1 function ip such that ip(x) $^(2/, x) 
for all x, with equality if and only if x = y, and 

Hi(x,DrP{y)) + (B{y)$. ti (y,y)). < 0. 

We define a function w € (Lip(T Ar )) m by setting 

Wi(-) = msoL{$i t i(y,-),ip + e}, Wj(-) = ®j,i(y, •) for j ^ i. 

Arguing as in the proof of Proposition 13.31 we see that it is possible to choose e > 
in a such a way that w is a critical subsolution. Moreover, since W{ agrees with ip + e 
in a neighborhood of y, there exist S > and an open neighborhood W of y such 
that u>i is of class C 1 in VF and 

Hi(x,Dwi(x)) + (i?(x)w(x)) i < —S for every x € W. 

Conversely, let assume (ii). According to Lemma |3.7| there exists w G %(0) such 
that Wi is smooth and strict in a neighborhood of y, in particular 

HifaDwiiyj) + [B(y)w(y)) i <0. 

In view of Proposition 13.81 we conclude that y A. □ 

Remark 3.10. Proposition 13.81 expresses the fact, roughly speaking, that the i— 
th component of a critical subsolution cannot be strict at y. However, since the 
supersolution test (|23p is void when D~~Ui(y) is empty, this fact cannot be directly 
used to prove the equivalence stated in Proposition 13.91 This is precisely the reason 
why we needed the regularization Lemma 13.71 

We proceed by proving a global version of the previous proposition. We give a 
definition first. 

Definition 3.11. Let v € H(0). We will say that v is strict at y if i>i is strict at y 
for every i £ {1, . . . , m}. We will say that v is strict in an open subset U of T N if 
it is strict at y for every y S U . 

Theorem 3.12. There exists v S %(0) which is strict in T N \A. In particular, the 
Aubry set A is closed and nonempty. 

Proof. Fix i € {1, . . . ,m}. We first construct a critical subsolution v* whose i-th 
component is strict in \ A. According to Proposition 13.91 for every y G \ A 
there exist an open neighborhood W y of y, a critical subsolution -w y and 5 y > such 
that 

Hi(x,Dw?(x)) + (B{x)w y (x)) t < -by for a.e. x G W y (24) 

The family {W y : y G \ ^4} is an open covering of \ A, from which we can 
extract a countable covering (W n ) n of \ A. For each n G N, let us denote by 
(w n ,<5 n ) the corresponding pair in H(0) x (0, +oo) that satisfies (f24"|) in W n . Up to 
subtracting to each critical subsolution w n a vector of the form k n 1 , we can moreover 
assume that 10^(0) = 0. Hence the functions w n are componentwise equi-Lipschitz 
and equi-bounded in view of Proposition 12.21 in particular the function 

v l (x) = Y, ¥ ™ n (x), xeT N 

n=l 
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is well denned and belongs to (Lip(T JV )) m . By convexity of the Hamiltonians, for 
almost every x G we get 

H>{x,Dvl(x)) + {B(x)v i (x)) i < + {B(x)w n (x))^ < 0. 

n=l 

Moreover, the above inequalities hold with —5k/2 k in place of almost everywhere 
in Wk, for every k G N. This shows that v* is a critical subsolution, strict in \A. 
Now set 

m 

v(x) = V — v*(x), xGT^. 

A similar argument shows that v is a critical subsolution that satisfies the assertion. 

If A = 0, by compactness we would have 7i(—5) ^ for some 5 > 0, contradicting 
the definition of the critical value c = 0. □ 

In view of Proposition 13.91 we have the following characterization: 

Theorem 3.13. Let y G T N . The following are equivalent facts: 

(i) y i A; 

(ii) there exists w G 7i(0) which is strict at y; 

(Hi) there exists w G H(0) and % G {1, . . . , m} such that Wi is strict at y. 

We end this section by extending to weakly coupled systems a result which is well 
known in the case of a single critical equation. 

Proposition 3.14. A = f] {y G T N : (<S(i)w)(y) = w(y) for every t > 0} . 

wew(o) 

Proof. Let us denote by A' the set appearing at the right-hand side of the above 
equality. Fix a point y G A and let w be any critical subsolution. For every fixed 
index i G {1, . . . ,m} : the function u* = &. : i(y, •) + Wi(y)t satisfies w ^ u* in 
and Wi{y) = Ui(y). Moreover, u* is a critical solution, hence it is a fixed point for 
the semigroup S(t) by Proposition 12.61 By monotonicity of the semigroup, we have 

w l (y) < (SCtJw)^) ^ (<S(£)u%(y) = u\{y) for every t > 0, 

hence all the inequalities must be equalities, in particular uS(i)w) (y) = Wi(y) for 
every t > 0. This being true for every i G {1, . . . ,m} and w G H(0), we conclude 
that y G A'. 

To prove the converse inclusion, we take a critical subsolution v which is strict 
outside A. According to Proposition IA.31 for every y A there exists t y > such 
that («S(tj,)v) (y) > v(y), that is y A'. □ 

4. Regularization 

The aim of this section is to show how a strict critical subsolution can be regular- 
ized outside the Aubry set. In the case of a single critical equation, it is known that 
such procedure can be performed in such a way that the output is a strict critical 
subsolution which is, in addition, of class C 1 on the whole torus, see [31 [I9l [20] . This 
result holds for Hamiltonians that are locally Lipschitz in (x,p) and strictly convex 
in p and the proof relies on the following two facts: first, any critical subsolution 
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is differentiable on the Aubry set and, second, its gradient is independent of the 
specific subsolution chosen. This latter rigidity property holds for weakly couples 
systems too, as we will show at the end of the current section. What prevents us 
to extend to systems the existence of C 1 strict critical subsolutions is the lack of 
information on differentiability properties of critical subsolutions on the Aubry set. 

We first deal with the regularization issue. The tools are not new and are mainly 
borrowed from \19\ I20|. However, we provide a proof for the reader's convenience. 

We start with a local regularization argument. 

Lemma 4.1. Let u G %(0) and assume that, for some r > ; 5 > and y G \.4 
and for every i G {1, . . . , m}, 

Hi(x, Dui(x)) + (-B(x)u(x)) . < —5 for a.e. x G £>2r(y)- 

Then, for every e > 0, there exists u £ G %(0) such that 

(i) ||u e - u||oo < e; 

(ii) u e = u in T N \ B r (y) ; 

(Hi) u e is of class C°° in B r j2{y) and satisfies 

Hi(x,Du £ i(x)) + (B(x)u £ {x)) i < -- 6 for every x G B r/2 {y). (25) 

Proof. Let 4> : — > [0, 1] be C°° function, compactly supported in B r (y) and such 
that <j) = 1 in B r / 2 {y)- Let (p n ) n be a sequence of standard moilifiers on WL N . For 
every n G N, we define a function w n G (Lip(T Ar )) m by setting 

wf(x) = (j)(x)(p n *Ui)(x) + (l — 4>{x))ui{x) for every x G and i G {1, . . . , m}. 

It is apparent by the definition that w n is of class C°° in B r / 2 (y) and agrees with u 
outside B r (y). Arguing as in the proof of Lemma 13.71 we see that it is possible to 
choose n large enough in such a way that w n is a critical subsolution and satisfies 
()25|) . Since w n ^ u in T^, the assertion follows by seting u e := w n for a sufficiently 
large n. □ 

We now prove the announced regularization result. 

Theorem 4.2. There exists a critical subsolution which is strict and C°° in T ra \^l. 
More precisely, for every critical subsolution v which is strict in T N \ A and for 
every e > 0, there exists v e G %(0) such that 

(i) ||v e - v||oo < e; 

(ii) v £ = v on A; 

(Hi) v e is C°° and strict in T n \ A. 

Moreover, the set of such smooth and strict subsolutions is dense in 7i(0). 

Proof. We first show how to regularize a subsolution which is strict outside the 
Aubry set. Let v be such a subsolution (given by Theorem 13. 12j) and fix e > 0. 
Since v is strict in \ A, there exists a continuous and non-negative function 
6 : -> R with 5" 1 ({0}) = A such that 

Hi(x, Dvi) + {B{x)\{x)) i ^ -S(x) in T N 

for every i G {1, . . . ,m}. Clearly, it is not restrictive to assume that the inequality 
5(x) < e holds for every x G T^. In view of Lemma |4.1| we can find a locally finite 
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covering (U n ) n of T N \A by open sets compactly contained in T N \A and a sequence 
(u n ) n of critical subsolutions such that each u n is C°° in U n and satisfies 

H^x.Du 1 }) + (B(x)u n (x)) i < -| 8(x) for every x G tf n , 

|u n (x) - v(x)| ^ <S(x) for every x G T^. (26) 

Set 

5 n := inf (5 (x) for every n G N 

and choose a sequence (r/ n ) n in (0, 1) such that, for every x G and n G N, the 
following holds: 

\H(x,p) - H(x,p')\ < for all p,p' G B K+ i with < rj n , (27) 

b 

where k denotes a common Lipschitz constant for the critical subsolutions, in par- 
ticular for all the u n . Last, take a smooth partition of unity ((fi n )n subordinate to 
(U n ) n and choose the functions u n in such a way that the quantities ||u n — v||oo, 
which can be be made as small as desired, satisfy 

Y ||u fc - v||oo H^fclL < Vn for every n G N. (28) 

ken 
u k nu n ^0 

That is always possible since the covering (U n ) n is locally finite. 
We now define v e : -> R m by setting 

oo 

v e (x) = Y^ i Pn(x)u n (x) in \ A and v £ (x) = v(x) on A. 

n=l 

By definition, v e satisfies assertion (ii) and is C°° in T N \ A. From (I26p we infer 
that |v e (x) — v(x)| ^ 5(x) in \ A, which shows at once that v e is continuous in 
and that it satisfies assertion (i). Moreover, by taking into account (|28p and the 
fact that Difk = 0, one obtains, for every x G U n and i G {1, . . . , m}, that 



fceN fceN 
u k nu n ^0 u k nu n ^0 



Dvf(x)- Y Mx)DuHx) =| Yl {^{x)-v{x))D Vk {x) < Vn , (29) 
l 

in particular 



\Dvf{x)\ ^ Vn + Y M*)\DuHx)\ < 1 + 



K. 



keN 
u k nu n ^0 

We infer that v e is Lipschitz-continuous in T^. In order to prove that v £ is a critical 
subsolution and is strict in \ A, it will be enough to show that 

HAx, Dvf(x)) + (B(x)v £ (x)) . ^ for a.e. x G T N , 

for alH G {1, . . . , m}. 

Since Dv £ (x) = Dv(x) for almost every x G A, it suffices to establish the claim 
in the complementary of A. To this aim, by recalling the definition of r\ n and by 
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making use of ([29]) and of Jensen inequality, we get that, for every x G U n and 
i G {l,...,m}, 

H^DvUx)) + (B{x)v £ (x)) i Hi(x, £ <p k (x)Dv£(xj) + -± 

fegN 

c/ fe ni7„^0 
+ £ ^(x)( J B(x)u fc (x)) j 

k&l 

u k nu n ^0 



u k nu n ^0 

<"3^) + -6<— ■ 

This concludes the proof of the first part of the statement. 

For the density, let u be any critical subsolution. Let v be a critical subsolution 
which is strict outside the Aubry set (whose existence is assured by Theorem 13. 12j) . 
Then, for any A G (0, 1), the function (1 — A)u + Av is a subsolution which is strict 
outside the Aubry set. This subsolution can therefore be regularized using the above 
procedure, giving a subsolution w which is strict and smooth outside the Aubry set. 
Moreover, both these steps can be done in such a way that || u — w II qq is as small as 
wanted. This establishes the density. □ 

We now additionally assume the Hamiltonians Hi to be strictly convex in p and 
derive some further information on the behavior of Clarke's generalized gradients of 
the critical subsolutions on the Aubry set. 

We start with a preliminary lemma. 

Lemma 4.3. Let y G A and let u 1 , -- , be critical subsolutions. Then, for all 
ie{l,... ,m}, 

i 

n d^Kx) + 0. 

k=l 

Moreover, it contains a vector pi which is extremal for all the sets d c u^(x) and which 
satisfies 

Hi{y,pi) + (B(y)u k (y)). = for every k G {1, . . . ,£}. 
1 t 

Proof. Let w = - ^ u fc G K(0) and let pi G d c Wi{y) be such that 

" k=l 

H i (y,p)+(B(y)w(y)).=0. 

Such a pi must exist because otherwise Wi would be strict at y. Note that, by strict 
convexity of Hi, the vector pi must be an extremal point of d c Wi(x), hence it is a 
reachable gradient of Wi . Let y n — >■ y be such that u\ is differentiable at y n for every 
k G {1, ...,£} and n G N, and 



1 1 

Dwi(y n ) = - ^2 Du i(yn) -> Pi- 



k=l 
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Up to extraction of a subsequence, we can assume that Du^(y n ) — > qk for all k G 
{1, ...,£}. Then one readily obtains, by Jensen's inequality, that 

1 £ 

= Hi(y, Pi ) + (S(y)w(y)). < j£ (fli(y,<fc) + (5(y)u fc (y)) .) < 0. 

fe=i 

Therefore, all the inequalities | {Hi{y, qk) + (B(y)u k (y))i) ^ summing to an equal- 
ity, we deduce, by strict convexity of Hi, that qi = ■ ■ ■ = qi = p- Moreover, since 

Hi(y, q k ) + (£(y)u fc (y)) . = for every ke {!,...,£}, 

and because of the strict convexity of Hi, one sees that pi is extremal, and thus 
reachable, for all the u\. □ 

We now extend the previous result as follows: 

Proposition 4.4. Let y G A. Then, for each i G {1, . . . ,m}, there exists a vector 
Pi G 1^ which is a reachable gradient of m at y for every u G H(0) and which 
satisfies 

Hi(y, P i) + (B(y)u(y)) i = 0. 

Proof. For each critical subsolution u, let us denote by the set of reachable 
gradients p of Ui at y that satisfy Hi{y,p) + [B(y)u(y)) i = 0. This set is not empty 
and compact. The proposition amounts to proving that 

n p f+ - 

ueH(o) 

If this were the case, by compactness we could extract a finite empty intersection. 
But this would violate the previous lemma. □ 



5. Rigidity of the Aubry set and Comparison Principle 

In this section we establish some interesting properties of the Aubry set and we 
provide uniqueness results for critical solutions. As we will see, such results will 
follow rather easily thanks to the information gathered so far. 

We start with a remarkable rigidity phenomenon that takes place on the Aubry 
set. 

Theorem 5.1. Let y G A and i G {1, . . . , m}. Then 

v(y) = $. ti (y,y) + Vi(y)l for every v G TL(0). 
In particular, v(y) — w(y) G Ml for any v, w G %(0). 

Proof. Take v G %(0) and set u := ^.^(y, •) +Vi(y)I. According to Proposition 
13.21 u is a critical solution satisfying v ^ u in T N and Vi(y) = Ui(y). By applying 
Proposition 12.31 with xq := y we get the assertion. □ 

Remark 5.2. On the other hand, the above property does not hold when y £ A. 
Indeed, the proof of Lemma 13. 71 shows that any critical subsolution v which is strict 
at y can be modified in such a way that the output is a critical subsolution all of 
whose components except one coincide at y with those of v. 
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Next, we derive a comparison principle for sub and supersolutions of the critical 
weakly coupled system (fT6|) which generalizes to our setting an analogous result 
established in [5] for Hamiltonians of a special Eikonal form, see Subsection 16.11 for 
more details. In particular, we obtain that A is a uniqueness set for the critical 
system. 

Theorem 5.3. Let v, u G (C(T N )) m be a sub and a supersolution of the critical 
weakly coupled system (|16p . respectively. Assume that 

m m 

Ai(x)vi(x) ^ Ai(x)ui(x) for every x G A, (30) 

8=1 i=l 

where A : A — > (M + ) m is a function satisfying 

m 

Aj(x) > for every x G A. 

i=i 

Then 

v(x) ^ u(x) for every x G T N . 
In particular, two critical solutions that coincide on A coincide on the whole T N . 

Remark 5.4. The above theorem also implies that two critical solutions u and 
v are actually the same if their i-th components coincide on A, for some fixed 
index i G {1, . . . ,m}. This is consistent with Theorem 15.14 which assures that this 
"boundary" condition amounts to requiring that u = v on A. 

Proof. In view of the density result stated in Theorem 14.21 the critical subsolution 
v can be approximated from below by a sequence of critical subsolutions that are, 
in addition, smooth and strict outside A. Indeed, just pick a sequence (w n ) n€ N such 
that ||w n — vjloo < n~ l and then define v n = w n — n _1 l which then verifies v n ^ v 
and ||v n — v||oo < 2n~ 1 . Clearly, each element of the sequence still satisfies the 
boundary condition ([30|h hence it is enough to prove the statement by additionally 
assuming v smooth and strict in \ A. 
Let us set 

M := max max (vi — Ui) 

and pick a point xq G where such a maximum is attained. By Proposition 12.31 we 
know that v(xo) = u(xo) + Ml. If xq A, then v\ would be a smooth subtangent 
to u\ at xo- The function u being a supersolution, we would have 

^ H^xq, D Vl (x )) + (B(x )u(x)) 1 = H^xq, D Vl {x )) + (B(x )v(x)) v 

in contrast with the fact that v is strict in T N \ A. Hence xo G A and by (|5U|) we 
get 

m m 

M^Ai(xo) = ^ Ai(x ) (vi(x ) - Ui(x )) ^ 0, 

i=l i=l 

i.e. M ^ as it was to be proved. □ 

Last, we show that the trace of any critical subsolution on the Aubry set can be 
extended to the whole torus in such a way that the output is a critical solution. 

Theorem 5.5. For any v G H(0), there exists a unique critical solution u such that 
u = v on A. 
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Proof. This assertion is derived as a consequence of Proposition 13.141 by setting 
u(x) = lim S(t)v(x) for every x G T^. 

t— >+oo 

Indeed, S(i)v = v on A for every t > 0. Since the functions {<S(i)v : t > 0} are 
equi-Lipschitz and non-decreasing with respect to t, we infer that S(t)v ^ u in 
T^. Last, u being a fixed point of the semigroup S(t), we get that u is a critical 
solution. □ 

6. Examples 

The critical value and the Aubry set for a weakly coupled system of the kind 
studied in this paper have, in general, no connections with those of each Hamiltonian, 
considered individually. This happens also in simple situations, see for instance 
Example 1.2 in |28| . In this section, we present a couple of examples where more 
explicit results may be obtained for the critical value and for the Aubry set. 

6.1. The setting of [5|. The first example we propose corresponds to the setting 
considered in [5]. Assume that all the Hamiltonians are of the form Hi(x,p) = 
Fi(x,p) — Vi(x), where: 

(a) Fi and Vi take non-negative values; 

(b) Fi is convex and coercive in p; 

(c) Fi(x, 0) = for all x G T N and i G {1, . . . , m}. 
Furthermore, assume that 

m 

n^- i ({o})^0. 
i=i 

Under these hypotheses, we claim that the critical value is (whatever the coupling 
matrix is) and that the Aubry set is nothing but 

TO 

•A=f)Vr\{0}). 

i=i 

Indeed, it is easily seen that the null function u° always belongs to H(0) under the 
first set of hypotheses. Therefore, 7i(0) and the critical value verifies c ^ 0. To 
see that there is actually equality, consider a point xq G nV^ _1 ({0}) and any (C 1 ) 
function u. An easy consequence of Proposition 11.41 yields that B(xq)u(xo) must 
have a non-negative entry, say i, hence 

Hi(x , Dui(x )) + (B(x )u(x )) i = Fi(x ,Dui(x Q )) + (B (x )u(x )) . ^ 0. 

Therefore, u cannot belong to a H(—e) for a positive e. The same argument can be 
adapted in the viscosity sense for any (non necessarily C 1 ) function. Therefore is 
the critical value. 

To prove that nUj -1 ({0}) is the Aubry set, first notice that, for every y 
PlT^~ ({0}), there exists an index j such that Vj(y) > 0. Then the j-th compo- 
nent of the null function u° is strict at y. In view of Theorem 13.131 we get the 
inclusion 

m 

AQP\Vr\{0}). 

i=i 

The opposite inclusion is obtained as previously. Take any u € H(0) and xo G 
n^ _1 ({0}). We will do as if u is differentiable at xq, but the argument carries on 

27 



in the general case using test functions and the viscosity subsolution property. At 
xq we must have 

(Fi(xo,Dui(x ))) + B{x )u(x ) < 01. 

V / i€{l,...,m} 

But this is only possible if Dui(xo) = for all i 6 {1, ... , m} and if B(xo)u(xq) = 0. 
Indeed, otherwise, B(xq)u(xo) will have a positive entry, which is impossible. In 
particular, the above inequality holds with an equality. Since this happens for any 
critical subsolution u, we get xq G A in view of Theorem 13.131 As a byproduct, 
this also establishes that at any point of A, any critical subsolution must take as 
value a vector belonging to 11. This is a particular case of Theorem 15.11 (and in 
this particular form, it also appears in a weaker form in [3]). 

6.2. Commuting Hamiltonians. In this second example we consider the case 
when the Hamiltonians are strictly convex and pairwise commute. If the Hamilto- 
nians are of class C 1 , that means 

for every i,j G {!■?••• >^}- If the Hamiltonians are only continuous, the commu- 
tation hypothesis must be expressed in terms of commutation of their Lax-Oleinik 
semigroup, see |llj for more details. We also make the additional assumption that, 
individually, all the Hamiltonians have as critical value. Then, we claim that is 
the critical value of the system as well (whatever the coupling is). 

Indeed, it is proved in [30] that the Hamiltonians have the same critical solu- 
tions. In particular, there exists a function u G Lip(T N ) satisfying 

H{(x, Du) = in for every i G {1, . . . , m} 

in the viscosity sense. Since the coupling is degenerate, we infer that the function 
u° = ut is a solution of 

Hi{x,Du°i) + (B(x)u°(x)) i = inT^ for every i G {1, m}. 

Therefore, the claim is a direct consequence of Proposition 12.111 Moreover, in this 
setting, we may localize the Aubry set of the system using those of the individual 
Hamiltonians. In order to do so, let us recall another result from [llj . 

Theorem 6.1. Let H±, ■ ■ ■ ,H m be pairwise commuting and strictly convex Hamil- 
tonians, with common critical value equal to 0. Then they have the same Aubry set 
A* . Moreover, there exists a common critical subsolution v which is smooth outside 
A* and strict for each Hamiltonian, i. e. 

Hi(x, Dv(x)) < for every x G \ A* and i G {1, . . . , m}. 

Using this theorem, we easily see that the inclusion A C A* holds. Indeed, the 
function ~v(x) := v(x)l is a critical subsolution for the system which is strict outside 
A*. 

We also note that, as in the previous example, u(y) £ Ml for every y G A and 
every u G H(0) in view of Theorem 15. 1[ 

A particular case of this example is when all the Hi are equal. In this case we get 
the more precise statement: 
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Proposition 6.2. Let H be a convex Hamiltonian and assume H\ = ■ ■ ■ = H m = H . 

Then A = A* . Moreover, all critical solutions of the system are of the form vl = u1 
where u is a critical solution of H . 

Proof. The inclusion A C A* can be proved arguing as above (note that we do not 
need the strict convexity assumption here). Let us prove the converse statement. 
Pick v G H(0) and set v := max(uj, i G {1, . . . ,m}). We claim that v is a critical 
subsolution for H. Indeed, let x G and p G D + v(x). Then v(x) = V{(x) for some 
i G {1, . . . ,m}. Since v ^ v j with equality at x, we get p G D + Vi(x). We now use 
the fact that v is a subsolution of the system to get 

H(x,p) ^ Hi(x,p) + (B(x)v(x)). ^ 0, (31) 

where the first inequality comes from the fact that 

m m 
(B(x)v(x)). = J2 b ij( X ) V j( X ) > ^2 b ij(x)vi(x) = 0, 

i=i j=i 

which holds true since bij{x) ^ and Vj(x) ^ Vi(x) for every j ^ i. Let us now 
assume that v is strict outside A. Then the right inequality in (|3ip is strict as soon 
as x ^ A, yielding that v is a subsolution for H which strict on the complementary 
of A. This proves that A* C A, hence A = A*. 

Let now u be a critical solution for the system. Then v := max(tij, i G {1, . . . , m}) 
is a critical subsolution for i?. Moreover, as 

u(x) = ui(x)t for every x £ A, 

we deduce that v = u\ on A Since A = A*, there exists a critical solution u for H 
such that {t = v on *4. Now the function u = ut is a critical solution of the weakly 
coupled system satisfying u = u on i. By the comparison principle, i.e. Theorem 
15.31 we conclude that u = u. □ 

Appendix A. 

In this appendix we want to give a proof of Proposition 12.71 

In what follows, a function u will be said to be semiconcave on an open subset U of 
either or x M + if, for every x G U, there exists a vector p x G M. N such that 

u(y) — u{x) ^ (p x , y — x) + \y — x\ u{\y — x\) for every y G U, 

where u is a modulus. The vectors p x satisfying such inequality are precisely the 
elements of D + u(x), which is thus always nonempty in U. Moreover, d c u{x) = 
D + u{x) for every x G U. By the upper semicontinuity of the map x t-> d c u(x) with 
respect to set inclusion, we get in particular that Du is continuous in its domain of 
definition, see [6]. 

Next, we prove the following 

Proposition A.l. Let T > and G : [0,T] x T N x R N ^ R be a locally Lipschitz 
Hamiltonian such that G(s, •) is a strictly convex Hamiltonian, for every fixed 
s G [0, T\. Let u(t,x) be a Lipschitz function in [0,T] x T N that solves the evolutive 
Hamilton- J acobi equation 

— + G(t,x,D x u) = in(0,T)xT N , (32) 
in the viscosity sense. Then 
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(i) for every < r < T , the function u is semiconcave in [r, T) x T^; 

(ii) if u(0, •) is semi-concave in T N , then the functions {u(t, ■) : t E [0,T)} are 
egm-semiconcaue. 

Proof. Since u is Lipschitz, up to modifying G outside [0, T] x x Br for a 
sufficiently large R > 0, we can assume that G is superlinear in p, uniformly with 
respect to (t,x). We are then in the setting considered by Cannarsa and Soner in 
[7] and item (i) follows from their results. 

Let us prove (ii). Let us denote by L(t, x, q) the the Lagrangian associated with G 
through the Fenchel transform and by uq the initial datum u(0, •). It is well known, 
see for instance [Bj, that the following representation formula holds: 

u(t,x) = inf (uo(£(0)) + [* L(s^(s),i(s))ds), (t, x) G (0, T) x T^, (33) 

where the infimum is taken by letting £ vary in the family of absolutely continuous 
curves from [0, t] to T N . Moreover, the minimum is attained by some curve 7, which 
is, in addition, Lipschitz continuous (actually, of class C 1 ), see [9]. 

We claim that there exists a constant k, only depending on G and on the Lipschitz 
constant of u in [0, T] x T N , such that H7H00 ^ «■ To this aim, we apply Proposition 
2.4 in [22] to the function u(t,x) and the curve s i-> (s,7(s)) to get 

^ w(s, 7(3)) =p t (s) + (pa.(s), 7(s)) for a.e. s G [0, t], (34) 

where m (pt{s),p x (s)) is a measurable and essentially bounded function on [0,t] 
such that 

(Pt(s),p x (s)) G 9 c n(s,^(s)) for a.e. s G [0,t]. 
By integrating (|34p and using the Fenchel inequality we get 

u(t,x) = u o {y(0)) + [ Pt(s) + {p x (s), A f(s))ds 

Jo 

< «o(7(0))+/ p t (a) + G(s,7(*),P iB (*))+^(a ) 7(a),7(*))da 

JO 

< «o(7(0)) + / L( S ,7( S ), 7 ( S ))ds, 

•/ 

where in the last inequality we used the fact that u is a (sub)-solution of the time 
dependent equation, i.e. 

p t + G(t,x,p) ^ for every {p t ,p x ) G d c u(t,x) and (t,x) G (0,T) x T^. 

Since 7 is minimizing, all the inequalities must be equalities, in particular we obtain 

j(s) G d p G(s,j(s),p x (8)) for a.e. s G [0,t\. (35) 

This proves the claim by choosing 

K:=sup{|g| : q G d p G(s, x,p), (s,x) G [0,7] xT^, |p| ^Lip(u;[0,T] xT w ) }, 

which is finite thanks to the convexity and the growth assumptions assumed on G 
with respect to p. 

Let us now fix t G (0, T), x±,X2 G T N , X G [0, 1] and set x = Xx\ + (1 — X)x2- Note 
that x\ = x + (1 — X)h and X2 = x — Xh for h = X\ — x<i- Let us denote by 7 a curve 
realizing the infimum in f|33f) for such a pair of (£, x), by K a Lipschitz constant for 
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L restricted to [0, T] x x B(0,2n) and by u a semi-concavity modulus for no- 
We get 

Xu(t, x\) + (1 — X)u(t, X2) — u(t, x) 



a(«o( 7 (0) + (1 - \)h) + J L(a, 7 (fl) + (1 - A)/i, 7 ( S ))ds) 
+ (1 - A) (n ( 7 (0) - A/i) + C L(s )7 (s) - A/i, 7 (s))ds) 



v no(7(0))+ ^ L( S , 7 ( S ), 7 ( S ))d Sy 
= An (7(0) + (1 - A)/*) + (1 - A)u (7(0) - Xh) - 1*0(7(0)) 
+ A(^ (l(s, 7 (a) + (1-A)fc, 7 (s)) - L( S)7 ( S ), 7 ( S )))d S 

+ (1 - A) jT (l( S , 7 (s) - Afe, 7 00) - L( s ,7( a ), 7 ( s )))d S 

< A(l - X)(lu(\xi - x 2 \) + tK\x\ - x 2 \), 
which proves the assertion. □ 

The result just proved will be applied to weakly coupled systems as follows: 

Proposition A.2. Let T > and u = (u 1 ,...,u m ) G (Lip([0,T] x J N )) m be 
a solution of the evolutionary weakly coupled system ([?]). Let i G {1, . . . ,m} and 
suppose that Hi is locally Lipschitz and strictly convex. Then, for all < r < T, 
the function m restricted to [t, T) x is semiconcave. Moreover, if, the initial 
condition Uj(0, •) is semiconcave, then the functions {ui(t, •) : t G [0, T] } are equi- 
semiconcave. 

Proof. The function m solves, for the given index i G {1, . . . , m}, a Hamilton-Jacobi 
equation of the kind ([32]) with 

G(t,x,p) = Hi(x,p) + (B(x)u(t,x)) v (t,x,p) G [0,T] xfx R N . 

The conclusion follows by applying Proposition IA.ll □ 

We are now ready to prove Proposition 12.71 

Proof of Proposition 12.71 We recall that, by convexity of the Hamiltonians, 
subsolutions to the critical system coincide with almost everywhere subsolutions. 
This fact will be repeatedly exploited along the proof. 

Assume first that t h-> S(t)u + tat is non-decreasing. Pick to > such that the 
map (t,x) 1 y S(t)u(x) is differentiable at (to,x) for almost every x G and 

dtS(to)u(x) ^ —a for a.e. x G T^. 

By the Lipschitz character of the map (t, x) 1— >• S(t)u(x) and Fubini's theorem, this 
holds true for almost every to > 0. Using the evolutionary equation, which is verified 
at every differentiability point of S(t)u(x), we deduce that, for every i G {1, . . . , m}, 

Hi^x,D(S(t )u).(x)^ + (B(x)S(t )u(x)) i < a for a.e. x G T N , 
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that is, S(to)u G H(a). This being true for almost every to > 0, the conclusion 
follows by stability of viscosity subsolutions. 

Let us now assume reciprocally that u G T-L(a). We first approximate each Hamil- 
tonian Hi with a sequence (H^)k of convex Hamiltonians that are, in addition, 
locally Lipschitz in (x,p) and strictly convex in p. This can be done by taking a 
sequence (pk)k °f standard mollifiers on 1^ and by setting 

|2 



fl?(x,p) = jf p fc (i/)fTi(i - y, P ) dy + (x,p) g T 



— S k (t)u n (x) ^ -(a + 3e)l for a.e. (i,sc) € (0,r) x T^, 



Note that, for each index i G {1, ... ,m}, H\ =4 if* in x as fc -> +oo. Let 
us denote by %k( a ) the set of a-subsolution of the weakly coupled system © with 
a = at and H\, . . . , iJ^ in place of Hi, ... , H m , and by S k the semigroup associated 
with the corresponding time-dependent equation ((7J). 

Next, we approximate u with a sequence of (u n ) n of functions that are component- 
wise semi-concave by setting 

ufix) = inf uAy) + n\y — x\ 2 for every x € and i = 1, • ■ • , m. 

Fix e > 0. A standard argument shows that, for n large enough, u n G 7i(a + e). 
Moreover, by the Lipschitz character of u n and by the local uniform convergence of 
(H*, . . . , to (H\, . . . , H m ), we also have that u n G Hk(a + 2e) for k sufficiently 
large. We now apply Proposition IA.2I to infer that the map (t, x) i-> S(t)u n (x) is 
semiconcave in [0, r] x for every r > 0. By using the fact that the gradient of a 
semiconcave function is continuous in its domain of definition and by choosing r > 
small enough, we get S(t)u n G T~Lk(a + 3e) for every ^ t ^ r. By exploiting this 
information in the evolutive weakly coupled system, we get 

d 
Wt 
i.e. 

S k (t + h)u n ^ Sfc(i)u n - h(a + 3e)l for every < t < t + h < r. 

Now, by the comparison principle for the evolution equation and by using the fact 
that the semigroup commutes with the addition of scalar multiples of the vector 1 , 
we obtain that t h-> Sk(t)u n — t(a + 3e)l is non decreasing. We now exploit the fact 
that 

S k (t)u n =t S{t)u n and S(t)u n =t S(i)u in M + x 

fc->+oo n—>+oo 

to infer that 

i i — ^ <S(i)u Tl — t(a + 3e)l is non-decreasing on [0, +oo). 

Being this true for every e > 0, we finally have that t h-> <S(t)u n — ial is non- 
decreasing on [0, +oo). 

The last assertion follows from the equivalence just proved, together with the fact 
that the semigroup S(t) is non-decreasing and commutes with addition of vectors 
of the form a 1 with a G R. □ 

The next result can be seen as a refinement of Proposition 12. 7\ where we assume 
inequality to hold only for one of the equations of the system. The proof follows via 
the same argument and is omitted. 
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Proposition A. 3. Let u G (Lip(T Ar )) m and assume that there exist an index i G 
{1, . . . , to}, an open subset U of T N and a € R st/c/i i/iai 

H^x, Dui(x)) + (i?(rc)u(x)) i ^ a /or a.e. x E U. 

Then for every open set V compactly contained in U there exists ry > suc/t i/iai 

(^S^u^x) ^ Ui(x) — at for every x € V and t € [0, ry]. 
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